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PREFACE 


The  in-flight  vibration  test,  carried  out  under  conditions  very  close  to  those  of  the 
future  operation  of  the  aircraft,  is  the  most  conclusive  test  to  ensure  the  safety  of  a  new 
type  of  aircraft,  since  it  does  not  involve  any  of  the  theoretical  assumptions  on  which 
predictive  calculations  arc  based.  While  flight  flutter  tests  arc  delicate,  its  application 
is  essential  if  theoretical  calculations  have  revealed  a  tendency  to  flutter  in  certain  areas 
of  the  aircraft,  even  if  a  considerable  safety  margin  is  predicted. 

Such  tests  can  be  carried  out  according  to  various  processes  which  arc  selected  in 
relation  to  the  size  of  the  aircraft  and  the  accuracy  aimed  at.  In  an  early  article  of 
the  Manual  on  Aeroclastieity  (Volume  IV.  Chapter  10,  1961),  the  methods  used  in  this 
area  about  twelve  years  ago  are  described  by  Messrs  M.O.W.Wolfe  and  W.T.Kirby.  In  a 
more  recent  supplement  (Revision  1969).  French  developments  relative  to  the  various 
means  of  excitation  and  the  various  methods  for  utilizing  the  data  collected  were 
reported  by  Mr  Piazzoli. 

This  is  a  new  addition  in  which  the  authors  review  the  testing  methods  used  in  the 
United  States.  This  document  will  undoubtedly  prove  extremely  useful  to  engineers 
and  specialists  desirous  of  being  informed  of  the  present  state-of-the-art  in  this  particular 
field  and  may  even  provide  guidelines  for  the  selections  which  they  will  be  called  upon 
to  make  in  the  future. 


R.MAZET 

Editor  of  the  Manual  on  Aeroclastieity 
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REUEM  DEVELOPMENTS /V  FLIGHT FLUTTER  TESTING  IN  THE  UNITED  STA  TES 


/.  INTRODUCTION. 

The  critical  development  schedules  and  high  performance  requirements  of  recently  developed  military  and  commercial  aircraft  in 
the  United  States  have  demanded  the  use  of  new  techniques  for  the  rapid  and  accurate  determination  of  flutter  characteristics  and  the 
establishment  of  safe  flight  envelopes.  Although  significant  advances  have  been  made  in  theoretical  flutter  analysis  methods  and  in  wind 
tunnel  testing  techniques  over  the  last  10  years,  it  has  been  only  in  the  recent  past  that  flight  flutter  testing  methods  have  shown 
significant  advances  in  the  state  of  the  art.  The'c  advances  arc  principally  the  result  of  more  sophisticated  analysis  methods  which  are 
now  possible  Because  large  high  speed  computers  have  lieen  dedicated  to  the  flight  test  operation  to  reduce  the  overall  flight 
demonstration  time  required.  Since  consideration  of  flutter  is  a  pacing  item  for  expanding  the  speed  envelope  for  any  new  airplane, 
considerable  effort  has  lieen  spent  to  reduce  the  number  of  flights  required  to  obtain  flutter  data. 

There  arc  several  different  flight  flutter  testing  techniques  currently  in  use  or  being  developed  in  the  United  States.  Most  of  the 
methods  employ  some  form  of  sinusoidal  excitation,  using  either  aerodynamic  vanes,  internal  mass  shakers  or  power  control  systems. 
The  more  promising  methods  all  seem  to  employ  fast  sweep  rates.  The  old  reliable  shake  and  stop  method  still  appears  to  be  in  use.  as 
docs  the  pilot  impulse  technique:  there  has  also  been  at  least  one  recent  use  of  ballistic  impulses.  The  use  of  random  excitation 
techniques  does  not  seem  to  have  gained  much  popularity  in  the  U.S.  as  yet:  however,  several  of  the  analysis  schemes  appear  to  have  the 
capability  of  handling  the  response  to  random  inputs. 

With  regard  to  data  reduction,  the  Kcnncdy-Paucu  method.  Reference  I.  or  modifications  of  it.  still  appears  to  have  merit  and  is 
being  used  by  several  companies.  Techniques  based  on  the  fast  Fourier  transform.  Reference  2.  have  also  gained  popularity  for  flight 
flutter  testing. 

This  paper  presents  a  few  comments  on  some  flight  flutter  testing  procedures  in  use  or  under  development,  as  reported  by  some  of 
the  major  U.S  aerospace  companies.  No  attempt  is  made  to  provide  a  comprehensive  review  of  flight  flutter  testing  philosophy  since  a 
great  many  papers  have  already  been  written  on  the  subject,  e.g..  References  3  through  7.  The  main  part  of  the  paper  is  devoted  to  the 
Grumman  model-matching  approach. 

NOT  A  TION 

bQ-  bs.  c0.  cs  Coefficients  used  in  the  representation  of  the  unsteady  generalised  aerodynamic  force 

fit!  Genera lired  force  due  to  a  vane  shaker  or  mass  shaker,  etc. 

fhIU  Approximate  generalized  force  obtained  by  applying  a  data  hold  to  a  sampled  force  signal 

dv  Displacement  at  the  shaker  location  in  the  mode  which  is  responding 

m  Generalized  mass  at  zero  airspeed  based  on  a  unit  modal  displacement  at  the  location  of  the  response  transducer 

n  Integer  used  to  identify  a  particular  sampled  value 

s  Laplace  transform  variable 

t  Time 

v  Flight  speed 

z  Z  transform  variable 


7 

ittl 

ijltl 

x,.\2 

XJ  =  Xr 


wn 

ftu 

fls) 

F(z) 


Coefficient  used  in  the  exponential  approximation  to  the  indicia!  aerodynamic  force 
Dirac  della  function 
Sampling  function 

Complex  conjugate  roots  which  determine  the  oscillatory  motion  characteristic  of  an  aeroelastic  mode 
Real  root  which  determines  the  exponential  motion  characteristic  of  an  aeroelastic  mode 
Sampling  interval 

Resonant  frequency  at  zero  airspeed 

Generalized  coordinate 

Laplace  transform  of  fill  l  | fl 1 1 1  «  f(s| 

Z  transform  of  ft  nrl  .Z'lftsll  *  Ftzl 


Z  transform  of  ft  nrl 


-1  SUMMARY  OF  FLIGHT  FLUTTER  TEST/SG  METHODS 


The  following  summary  describes  the  most  recent  methods  used  by  some  of  the  major  United  States  aerospace  companies 
McDonnell  Mrcraft  Company 

Excitation  is  accomplished  through  the  stabilator  and  aileron  actuators  by  mean;  of  an  electrical  signal  which  feeds  directly  into 
the  servo  of  the  power  cylinders.  The  exciter  can  l>e  set  either  on  automatic  linear  frequency  sweep  or  be  manually  controlled  by  the 
pilot.  In  the  latter  case,  constant  frequency  shutoff  operation  of  the  exciter  may  be  employed  to  obtain  decay  records.  A  fully 
automated  analysis  for  determining  frequency  and  damping  from  a  transfer  function  has  been  developed  based  on  a  modification  to  the 
Ken.tedy-Pancu  method.  These  data  are  stored  in  a  computer  as  a  function  of  altitude  and  Mach  number.  Tracking  and  curve  fitting  of 
frequency  and  damping  as  a  function  of  dynamic  pressure  are  employed  to  identify  the  critical  mode;.  At  each  speed  after  the  first  two, 
th-  flutter  margin.  Reference  S.  is  calculated  at  constant  altitude  or  constant  Mach  number  as  desired.  A  cathode  ray  display  of  any  of 
the  derived  quantities  is  available  on  line  in  ewcntially  real  time.  The  accuracv  of  the  determined  frequency,  damping,  and  flutter 
margin  has  been  checked  using  theoretical  and  wind  tunnel  model  data. 

Douglas  Aircraft  Company. 

The  DC-10  airplane  has  aerodynamic  vanes  installed  on  all  surfaces  for  flutter  excitation.  The  response  transducer  signals  are 
digitized  and  stored  on  magnetic  tapes  in  the  flight  test  aircraft,  and  are  abo  telemetered  to  the  ground  facility.  The  ground  facility 
provklcs  instant  replay  of  the  last  45  seconds  of  data  on  live  analog  display  and  also  provides  for  post  flight  digital  analysis.  One  notable 
method  of  post  flight  analysis  involves  a  least  squares  fit  of  the  aircraft  transient  response  to  pitot  inputs.  The  analysis  follows  the 
technique  given  in  Reference  9  except  that  a  two  step  iteration  procedure  is  used.  The  zero  offset  amplitude,  and  phase  angle  art  found 
in  the  first  step  by  a  direct  least  squares  analysis:  a  Taylor  expansion  and  ieast  squares  fit  b  used  to  find  the  damping  and  frequencies  in 
the  second  step.  This  two  step  procedure  is  repeated  several  times,  converging  on  all  terms  to  completely  define  the  amplitude,  phase 
angle,  frequency,  and  damping  of  each  mode  in  the  response.  The  operator  has  a  choke  of  searching  for  one.  two.  three  or  four  degrees 
of  freedom 

It  has  been  found  that  reruns  with  minor  changes  in  the  start  and  stop  times  give  repetitive  results  for  "true  data"  but  results 
derived  from  noise  are  different  for  each  modified  run.  The  computer  b  programmed  to  ignore  any  amplitudes  that  do  not  exceed  a 
certain  percentage  of  the  maximum,  thus  preventing  it  from  trying  to  converge  on  very  low  level  signals. 

t.nckheeJ-Cahfornia  Company 

The  Lockheed  technique  utilizes  the  response  to  variable  frequency  sinusoidal  excitations  to  evaluate  the  dynamic  acroelastic 
stability  (modal  damping)  of  the  aircraft  structure.  The  method  b  based  upon  the  theoretical  response  of  a  single-degree-of-freedom 
spring-mass-damper  system  when  excited  by  a  sinusoidal  force  of  constant  amplitude,  but  of  linearly  increasing  or  decreasing  period. 
The  method  was  first  developed  by  l)r.  E.  A.  Bartsch  of  Lociheed-Califomia  Company  in  conjunction  with  flutter  test*  on  the  F-104 
fighter  and  the  Electra  transport  aircraft  beginning  about  1957.  The  distinguishing  feature  of  the  method  b  its  reliance  on  the  frequency 
sweep-rate  effect.  The  sweep  effect  produces  a  shift  in  the  response  frequency  at  resonance  which  increases  with  increasing  sweep  rate 
and  with  decreasing  damping.  By  evaluating  the  shift  in  resonant  frequency  between  an  up-sweep  and  a  down-sweep,  relative  to  (he 
average  response  frequency,  a  very  sensitive  measure  of  damping  b  obtained  for  all  critical  modes. 

in  practice,  the  computerized  version  of  the  analvsb  technique  uses  flight  recorded  structural  responses  to  symmetric  and 
antbymmctric  variable  frequency  excitation  forces.  These  forces  are  generated  aerodynamically  on  the  L-lOi  I  wing  by  wing  tip  vanes 
and  on  the  horizontal  tail  by  two  tail  tip  vanes.  The  analysb  data  consbt  of  time  series  measurements  of  the  appropriate  pair  of  vane 
excitation  forces  together  with  10  to  20  wing  ur  tail  response  measurements.  The  program  tabulates  the  time  series  of  the  excitation 
forces  and  the  responses,  as  well  as  the  oulput/input  ratios  and  the  instantaneous  frequency  sweep  rate.  Rots  of  these  quantities  versus 
time  arc  provided  together  with  plots  of  the  output/input  ratios  versus  frequency.  Tlie  peak  response  frequencies,  frequency  shifts,  and 
amplitudes  are  then  evaluated  to  determine  the  frequency  and  damping  of  the  critical  acroelastic  modes  for  all  flight  conditions  of 
interest. 

General  Dynamics  Corporation: 

Flight  flutter  testing  follows  the  practice  of  detecting  major  response  frequencies  and  obtaining  the  damping  of  these  modes  at 
successive  speed  increments.  The  variation  of  damping  with  speed  is  monitored  via  telemetry  to  allow  safe  expansion  of  the  allowable 
flight  envelope. 

Excitation  b  provided  by  oscillating  the  control  surfaces  if  the  frequency  response  of  the  actuating  system  b  satisfactory  over  the 
frequency  range  of  interest.  Otherwise,  excitation  b  supplied  by  either  aerodynamic  «abs  or  oscillating  mass  shakers.  The  tabs  and 
shakers  are  driven  hydraulically  with  provisions  for  both  amplitude  and  frequency  control.  Automatic  frequency  sweep  capability  with 
sweep  rale  proportional  to  frequency  is  utilized.  A  typical  time  for  frequency  sweep  from  0  to  40  Hz  b  one  minute. 

Accelerometers,  strain  gages,  and  position  pickups  are  used  to  measure  the  airframe  response  to  excitation  The  output  of  selected 
pickups  b  plotted  on  x-v  plots  during  the  frequency  sweeps.  Major  resonances  observed  on  these  plots  ate  noted  and  corresponding 
damping  values  obtained  from  skakt-snd-stop  tests. 

The  Hoeing  Company 

The  following  paragraphs  describe  the  method  which  Boeing  proposed  for  use  on  its  SST  airplane. 

Aerodynamic  vanes  located  at  the  wint  tips  of  the  SST  were  selected  for  in-flight  excitation  of  wing  modes  which  were  of  primary 
concern  from  a  flutter  standpoint.  The  wing  tip  aerodynamic  vane  «i>  chosen  because  it  provides  adequate  force  for  the  least  weight. 
The  choice  was  based  on  the  results  of  analytical  studies  in  the  subcri'ical  speed  range  indicating  the  wing  tip  to  be  an  effective  location 
from  which  to  excite  the  wing  flutter  mode.  Maximum  frequency  for  the  aerodynamic  vane  shaker  was  to  be  approximately  20  Hz 


Response  of  the  airplane  to  a  rapid  sinusoidal  sweep  was  to  lie  measured  by  accelerometers  installed  in  the  wing  and  body.  The 
outputs  were  to  l>c  recorded  on  magnetic  tapes  and  digitized  Indore  undergoing  a  fast  Fourier  analysis.  The  output  was  to  appear  as  a 
frequency  spectrum  plot,  for  each  sweep  and  transducer,  obtained  from  the  ratio  of  Fourier  transforms  of  output  to  input.  These  data 
allow  calculation  of  the  modal  damping  for  each  sweep  by  the  Kennedy-Pancu  method. 

THE  MODE!.  MA  TI  WNG  METHOD 

In  1 968.  Grumman  decided  to  purchase  a  large  computer  facility  for  use  exclusively  by  the  flight  test  organization.  The  purpose  of 
this  investment  was  to  effect  a  major  reduction  in  the  time  taken  to  flight  test  a  new  airplane.  The  reduction  was  to  be  accomplished  by 
shortening  all  phases  of  the  lest  program:  one  of  the  most  significant  of  these  was  the  linn-  spent  at  each  test  point  during  the  high 
speed  build-up  to  acquire  nutter  data. 

When  the  decision  was  made  to  use  a  computer,  v/e  had  the  opportunity  to  develop  a  new  technique  tliat  could  provide  frequency 
and  damping  data  in  a  relatively  short  time.  A  review  of  Grumman  methods  in  use  at  that  time.  Reference  10.  showed  that  a  dedicated 
computer  could  readily  reduce  data  reduction  time:  however,  data  acquisition  time  could  not  lie  significantly  altered.  A  technique 
which  reduces  data  analysis  time  and  is  compatible  with  relatively  fast  data  acquisition  is  dcscrilied  in  Reference  1 1  and  became  the 
basis  of  the  current  Grumman  procedure.  The  initial  implementation  of  this  technique  was  accomplished  by  Astrodata.  Inc.  under  the 
direction  of  R.  G.  Madsen. 

The  technique  can  best  lie  described  as  "model  matching".  An  analysis  model  of  known  functional  lorm  represented  by  a 
difference  equation  is  programmed  in  the  computer:  constants  are  de  iieil  which  cause  the  model  response  to  lie  identical,  within  some 
prcscrilicd  degree  of  accuracy,  to  the  aircraft  response  to  an  input  force:  these  constants  arc  then  used  to  determine  the  damping  and 
frequency  of  the  models).  The  matching  is  done  for  the  aircraft  response  data  over  a  limited  frequency  range.  Each  significant 
frequency  hand  of  data  is  anal)  red  separately  and  a  different  model  match  is  computed. 

While  different  types  of  force  inputs  can  lie  used  with  this  method,  it  is  necessary  that  the  force  by  such  that  it  causes  the  modes 
in  question  to  lie  excited  to  reasonable  signabnoise  ratios  and  that  it  lie  related  tr  some  transducer  signal  by  a  transfer  function  of 
known  form.  In  order  to  keep  the  time  for  data  acquisition  to  a  minimum,  a  force  is  chosen  which  has  a  frequency  cuntent  consistent 
with  the  modes  of  interest  and  has  a  duration  of  application  only  long  enough  to  give  the  modes  timr  to  respond  to  practical  levels. 
Grumman  has  chosen  a  relatively  fast  frequency  sweep  for  its  method.  Sweeps  from  S  to  50  Hz  in  1 5  seconds  arc  typical.  This  rate  is 
fast  enough  to  get  the  test  data  quickly  hut  slow  enough  to  insure  that  the  analog  records  still  contain  the  familiar  resonance  peaks  of  a 
steady  state  response  plot.  This  latter  feature  is  desirable  to  permit  the  test  engineers  to  gain  additional  understanding  of  the  response 
characteristics  by  studying  the  analog  traces. 

The  derivation  of  the  difference  equations  for  the  analysis  model  will  now  lie  presenter!.  This  will  lie  followed  by  a  description  of 
the  associated  computer  program  and  the  methods  employed  to  gain  confidence  in  the  approach.  Finally,  some  results  from  a  recent 
flight  test  program  will  lie  shown 

J.  /  Derivation  of  the  Difference  Equation 

For  simplicity .  the  derivation  of  the  required  difference  equation  will  lie  given  for  a  single  mode  model.  The  comparable  derivation 
for  higher  order  models  would  only  require  inure  elaborate  algebraic  procedures 

We  begin  with  the  generalized  equation  of  motion  for  the  mode. 


i  rcj^f  + J  V  |1S£+  !v2c^+-i  -if (l|!+ v:c  £tl|)|e">"  ‘''dt,  =£lflt 
"  m  in  m  i  '  m 


The  first  two  terms  on  the  left  are  inertial  and  structural:  the  last  three  are  aerodynamic.  The,  form  of  the  unsteady  aerodynamics  has 
been  chosen  because  i»  is  applicable  to  nanharnomic  motion.  A  single  xponcnlial  term  lias  been  used  in  the  approximation  of  the 
indicia!  aerodynamic  response  for  simplicity  A  description  of  the  derivation  of  equations  of  this  type  for  a  general  system  of  modes  Is 
given  ir.  Reference  1 2. 

Equation  ( 1 1  describes  the  continuous  response  of  the  inode  starling  from  rest.  7  o  obtain  the  transfer  function  between  input  force  and 
acceleration  response,  the  Laplace  transform  is  taken 

a-  •  (snis2  - 

'*  1  in  I s-A  j  )(s-kslts-X  j)  0vfi'> 

where  |  and  \i  arc  the  complex  conjugate  roots  associated  with  the  frequen  y  and  damping  and  is  a  real  root  associated  with  the 
judicial  aerodynamics  It  can  be  shown  that  the  Vs  arc  roots  of  the  cubic  poly  nominal 
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The  olijcclivc  of  the  model-matching  program  is  to  solve  for  Hie  roots  of  Equation  (2)  using  a  difference  equation  which  relates 
sampled  values  of  acceleration  and  force  The  sampled  values  are  obtained  by  feeding  analog  response  signals,  telemetered  from  the 
aircraft,  to  an  analog  to  digital  converter  in  the  ground  station  and  then  passing  this  digital  data  directly  to  the  computer. 

The  differemc  equation  is  derived  In  using  the  theory  of  santplcd-tlaia  svstems.  Reference  Id  The  functional  representation  of 
the  sampling  process  is 


5  j  r  1 1  =  L  At  1  nn 


The  sampled  force  is  then  simply  f(t)  5y(t). 

The  Laplace  transform  of  this  sampled  force  is 

oo 

2  f(nr)  e'nn  . 

n=0 

Now,  a  hold  or  data  reconstruction  function  is  required  to  convert  this  sampled  function  back  into  one  wliich  is  essentially  the 
same  as  the  original  continuous  function  f(t).  Since  the  conversion  is  not  exact,  we  will  call  the  reconstructed  signal  fj,(t).  The  hold 
function  used  in  the  Grumman  model  is  called  a  polygonal  hold;  its  transfet  function  is  given  by: 

e^H-e'")2 

rs2 

We  illustrate  the  overall  process  in  Figure  1(A).  Using  Z  transforms,  the  desired  difference  equation  can  be  derived  for  the 
sampled-data  model  shown  in  this  figure.  Letting  z  =  e75, 

OO 

2  f(nr)z'n  =  F(z) 

n=0 


2  £(nr)z'n  =  Z  (z) 
n=0 

it  follows  that 

F(z)T,(z)  =  2(z) 


rs2  (s-Xj  )(s-X2){s-X3)_ 

If  the  indicated  Z  transform  is  carried  out,  we  hare 


a4  +  (as  -2a4)  z'1  +  (as  -2a4)  z'2  *  a;  z'3 


I-alZ*‘-*2*'2-»3Z*3 


a,  =  eXlT  +  eX2r  +  cX3r 

,2  =  -e<xl +  X2>r-«,X1 +x3)T-e,x2  +  X3,r 

a3  =  e(Xl +X2  +  X3)r 

*•  *  £  a,.,. [>*' '  «l  *  ■  ^3'  *  l»2  *  «>.3  ■  ‘  1 1  *  »2>] 

•!  ■  b  W  "y ’  Tl»l  -  ■ h’ *  'a.  *  a2  *  ,»J .  1, )  .  #2 'h>'  .  ,kx2  ■  x,,j 

By  taking  the  inverse  Z  transform  of  Equation  (4),  the  difference  equation,  in  the  sampled  time  domain,  relating  past  and  present 
values  of  the  system  input  and  acceleration  response  is  obtained: 

a4f(nr)  v  (as-2a4)  f(n-l  )r  +  (a4-2as)f  (n-2)r  +  ajf(n-3)r  =  £(nr)  -«,  {(n-1  )r  -a2  {<n-2)r  -a3  {(n-3)r.  n  =  3. 4.  S . 

This  equation  can  be  written  in  matrix  form  as  follows: 

{(2r)  i(r)  &0)  |f(3r)-2f(2rt+fir)|  (f(2r)-2f{r)+fl0)|  a,  t(3r» 

«3r)  j(2r)  *(r)  |f(4r)-2f(3rH-fi2r)|  |f(3r>-2f(2r>+f<r)|  a2  l<4r) 


j(i-l)r  j(i-2)r  £<i-3)r  |f(ir)-2f(i-l)r+f(i-2)r]  (fli-l)r-2f| 


i-2)r+f(i-3)rj 
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This  equation  can  be  solved  by  least  squares  for  the  unknown  coefficients  a  |  through  a$  by  letting  i  take  on  a  sufficient  number  of 
values  over  the  dat  .gnient  of  interest.  The  roots  Xj.  Xa.  and  X^.  can  be  obtained  from  aj.  a-,,  and  a^:  the  remaining  unknowns. 
34  and  35,  can  be  ustJ  to  obtain  X  and  dv/rm. 

The  derivation,  to  this  point,  has  assumed  that  f(t)  is  such  that  a  transducer  can  measure  the  applied  force  directly.  Actually,  the 
shakers  used  in  flight  flutter  testing  typically  do  not  satisfy  this  assumption.  To  avoid  this  difficulty.  3  transfer  function  relating 
transducer  signal  to  force  is  assumed  and  the  difference  equations  are  derives)  using  a  model  similar  to  Figure  1(B).  For  example,  a  vane 
shaker  has  inertial  and  unsteady  aerodynamic  force  terms: 

fv(t)  =  m%j3  +  vbv^  +  v-cv^  +  Aj  Jvb^  (tj )  +  v‘cv^fl  (t,  )J  eVMl '  dt  j 

The  symbology  in  this  equation  corresponds  to  that  in  Equation  ( 1 )  but  a  “v"  subscript  refers  to  the  vane  and  p  refers  to  vane  angular 
position.  The  transfer  function  between  fv(s)  and  (3(s)  is: 


fv(s)  =  rnv< s  -  f  |  )(s-  falls  -  $3) 


1 

I5  +  7y) 


P(  s) 


where  the  f's  are  the  roots  of  a  cubic  polynomial  simitar  to  (3).  If  the  aerodynamics  are  quai-steady  instead  of  unsteady,  as  would  lie 
the  case  when  the  vane  is  very  small  rebtise  to  the  surface,  we  would  write 

fv  (s)  =  m  (s  -  f.  )(s  -  {■■>)  i5s) 
qs  v  1  _ 

If  the  shaker  is  inertbl  only,  with  no  aerodynamics,  we  would  write 


fv  (s)  = 


A  difference  equation  would  now  be  derived  which  would  not  only  have  the  unknowns  aj  through  aj  as  in  Equation  (S).  but  would 
contain  additional  unknowns  corresponding  to  the  f's  and  yv. 


3. 2  Description  of  the  Computer  Program  and  Automated  Telemetry  Station 


The  Grumman  computer  program  for  model  matching  is  written  with  the  assumption  that  the  forcing  function  will  be  a  sinusoidal 
frequency  sweep  from  low  frequency  to  high.  The  program  takes  the  sampled  shaker  position  signal  and  continuously  calculates 
frequency.  The  test  engineer  provides  a  starting  and  stopping  frequency  for  each  data  segment  to  be  analyzed.  The  computer  program 
finds  the  data  sequence  within  this  segment,  fitters  it.  and  stores  the  information  into  the  matrix  array  described  in  the  bst  section. 
Before  storing  the  data  into  this  matrix,  a  decimation  is  accomplished  which  reduces  the  number  of  data  points/cycle  to  five  at  the  high 
frequency  end  of  the  band  being  analyzed.  This  reduction  is  done  to  reduce  computation  time. 

Options  are  provided  to  use  any  of  three  types  of  digital  filters  with  various  corner  frequencies  and  roll-off  characteristics.  Options 
are  also  provided  to  use  matrices  conforming  to  2nd.  3rd.  4th.  5th.  and  6th  order  models.  The  program  then  solves  the  matrix  equation 
by  inversion,  solves  for  the  frequency  and  damping  using  a  polynomial  solution  method,  and  displays  the  results  on  a  cathode  ray  tube 
together  with  appropriate  identification  data. 

The  program  provides  for  five  separate  frequency  sweeps  at  each  test  point.  These  correspond,  typically,  to  symmetric  and 
antisymmetric  wing  and  horizontal  tail,  and  antisymmetric  fin.  For  each  sweep,  provision  has  been  made  to  obtain  solutions  for  as 
many  as  seven  record  segments  with  separate  model  orders  for  each.  Any  of  approximately  20  response  channels  can  be  used  to  supply 
the  data  for  the  35  solutions.  A  primary  and  secondary  transducer  are  used  for  each  solution.  Normally,  the  test  engineers  supply  the 
necessary  clues  to  the  computer  in  advance  of  the  test  point  by  extrapolating  from  previous  data.  If.  after  the  sweeps  are  analyzed,  it  is 
found  that  some  of  the  clues  are  not  correct,  changes  can  be  made  while  the  flight  is  still  in  progress 

The  data  can  be  analyzed  in  nearly  real  time*  in  one  mode  of  operation.  This  mode,  however,  requires  3  large  amount  of  computer 
storage  and  prohibits  other  disciplines  (stability  and  control,  engine  performance,  etc.)  from  obtaining  their  test  data  in  real  time.  The 
preferred  mode,  for  high  speed  test  points,  is  one  wherein  the  five  shaker  sweeps  are  done  quickly,  one  right  after  the  other,  while  the 
data  is  digitized  and  stored  directly  on  a  disk.  After  all  test  objectives  are  completed,  the  aircraft  slews  down  and  waits.  The  flutter  disk 
is  then  called  up  and  the  data  analyzed  in  less  than  five  minutes. 

The  data  flow  through  the  Automated  Telemetry  Station  (ATS)  for  flutter  testing  is  shown  in  Figure  2.  In  the  ttr 
DEMODULATOR,  the  data  from  the  microwave  stream  are  demodubted  into  26.500  words  per  second  of  pulse  code  (PCM)  data  and 
two  carrier  bands  each  carrying  1 4  subcarriers  of  FM  flutter  test  information.  These  three  data  streams  arc  recorded  on  magnetic  tape  to 
form  a  complete  record  of  all  raw  telemetered  data  acquired  during  a  flight. 

The  PCM  signal  is  aligned  in  the  BIT  SYNCHRONIZER,  and  transferred  to  the  PREPROCESSOR.  The  28  FM  subcarriers  are 
further  demodulated  in  DISCRIMINATORS  to  obtain  the  raw  analog  flutter  data. 

The  1 2-bit  digital  output  words  from  the  Analog  -  Digital  Converter  (ADC)  are  passed  to  the  PREPROCESSOR  in  parallel  with  the 
PCM  data.  The  PREPROCESSOR  performs  the  following  functions: 


•  Converts  PCM  data  to  parallel  format  from  its  serial  telemetry  format 

•  Applies  calibration  to  the  data 

•  Buffers  the  data  for  transfer  to  the  Computer  (CPU)  at  0.1  second  intervals 


NHfly  rwl  time  me#n*.  to*  ***mpfe.  t*tt  the  tr»  to  for  Xv?n  record  ttgrrxnti,  pnmtry  *nd  *tcond**Y  *  *n#dwc?<t.  wnll  be  diKJftyed  th*n  30  tec.  after  the  completion  of 
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Records  alt  data,  after  calibration,  on  magnetic  tape  in  a  format  comparable  with  the  CPU 
Transfers  to  the  CPU  up  to  1 5,000  words  per  second  of  PCM  and  digital  data. 


• 

The  computer  is  a  CDC  6400  with  98.000  60-bit  words  of  central  memory  and  a  disk  storage  unit  capable  of  storing  16,500,000 
60-bit  words.  The  command  logic  of  this  computer  enables  it  to  serve  three  Data  Analysis  Station  (DAS)  terminals  simultaneously. 
During  flight  testing,  one  (or  more)  of  these  DAS  terminals  is  dedicated  to  the  data  stream  from  each  aircraft.  These  terminals  are 
shared  in  turn  by  various  test  disciplines,  so  that  more  than  one  aircraft  can  fly  simultaneously,  and  each  aircraft  may  acquire  data  to 
satisfy  various  test  objectives  on  the  same  flight. 

In  parallel  with  the  digital  data  flow,  16  of  the  FM  discriminator  outputs  are  passed  to  two  8-channel  BRUSH  RECORDERS. 
These  displays  arc  monitored  for  safety  of  flight  during  acceleration  to  the  next  eneviope  expansion  point  and  to  evaluate  data  quality 
which  may  influence  program  results  (e.g.,  telemetry  interruption  during  a  sweep). 

3. 3  Results  of  Model  Matching  When  Applied  To  Theoretical  Response  Data 

Prior  to  using  the  model  matching  program  on  a  new  airplane,  a  rather  elaborate  checkout  was  undertaken.  The  program  was  first 
used  to  analyze  digitally  generated  response  data  for  a  spring-mass-damper  system  with  six  degrees  of  freedom  and  then  it  was  used  to 
analyze  digitally  generated  response  data  for  a  two  mode  aeroelastic  system. 

In  the  analysis  of  the  six  degree  of  freedom  system  it  was  found  that  for  a  noise  free  sign-l.  the  model  matching  technique  would 
give  results  which  were  almost  exact  for  one  or  two  modes.  When  white  noise  was  introduced  in  the  response  signal  (noise  meaning 
response  which  was  not  due  to  the  known  force  input),  the  accuracy  was  somewhat  degraded.  It  became  apparent  during  this  early 
checkout  that  in  a  noise  envirionment  it  was  desirable  to  keep  the  record  segment  as  short  as  possible  and  confined  to  the  area  where 
the  response  signal  reached  a  peak  value. 

The  analysis  of  the  two  mode  aeroelastic  system  was  undertaken  to  make  a  quantitative  evaluation  of  the  ability  of  the  model 
matching  program  to  obtain  frequency  and  damping  values  from  data  generated  during  a  fast  sinusoidal  sweep.  The  influence  of 
unsteady  aerodynamics,  which  manifests  itself  in  the  modal  response  and  in  the  generation  of  the  forcing  function,  was  included  in  the 
evaluation  as  well  as  the  effect  of  atmospheric  turbulence 

The  theoretical  model  chosen  was  one  with  a  9.2  Hz  bending  mode  and  a  28.9  Hz  torsion  mode.  The  basic  parameters  were  chosen 
so  that  an  explosive  flutter  condition  occured  at  1065  kts.  The  unsteady  aerodynamic  forces  were  represented  using  an  exponential 
approximation  to  Wagner’s  indicia!  lift  function.  One  term  was  used  in  the  approximation:  6<s)  -  l-.Se'-‘34s.  A  vane  shaker  was  located 
at  the  trailing  edge  of  the  surface  and  its  hinge  line  unbalance  was  made  variable  so  that  the  frequency  at  which  the  aerodynamic  and 
inertial  forces  canceled  each  other  (the  crossover  frequency)  could  be  varied.  The  same  approximation  to  Wagner’s  indicia!  lift  function 
was  used  for  the  vane  aerodynamics. 

A  quasi-steady  gust  loading  was  applied  to  the  model  to  simulate  atmospheric  turbulence  in  the  following  manner:  a  computer 
program  which  generates  random  white  noise  was  used  as  a  downwash  input  to  a  first  order  differential  equation  which  shaped  this 
downwash  into  a  reasonable  approximation  to  the  gust  spectrum  shown  in  Reference  14.  This  gust  was  then  multiplied  by  the 
appropriate  constant  which  converted  it  to  a  force  which  was  apptied  at  the  20%  chord  of  the  surface.  The  rms  value  of  gust  intensity 
was  treated  as  a  variable. 

The  damping  and  frequency  curves  for  the  system  are  shown  in  Figure  3.  They  were  obtained  by  taking  the  Laplace  transform  of 
the  homogeneous  equations  and  solving  the  resulting  characteristic  polynomial  for  the  roots.  Forced  response  curves  for  the  system  are 
shown  in  Figures  4-10.  These  were  calculated  by  numerically  integrating  the  equations  of  motion  treating  vane  position  as  the  following 
known  input  function:  0  =  I-cos<  150e-^*).  The  velocity  response  at  a  specific  point  on  the  surface  and  the  vane  position  were  written 
on  a  magnetic  tape  using  the  same  time  increment  as  that  used  by  the  analog/digital  converter  for  telemetered  signals.  This  tape  was 
then  analyzeo  with  the  model  matching  program. 

Frequency  and  damping  values  obtained  with  the  model  matching  program  are  summarized  in  the  tablets  included  with  Figures  4 
through  10.  Three  air  speeds  are  presented:  I)  400  kts.  which  is  a  relatively  low  speed  for  the  system  and  one  wherein  very  little 
aerodynamic  coupling  exists  between  modes;  2)  840  kts.  which  is  sufficiently  high  to  insure  significant  aerodynamic  coupling  while 
being  well  short  of  the  coalescence  speed  and  flutter  speed;  and  3)  1000  kts.  which  is  the  coalescence  speed  and  is  just  below  the  flutter 
speed.  At  each  of  these,  results  are  presented  for  noise  free  data  and  for  data  which  includes  atmospheric  turbulence.  The  turbulence 
levels  chosen  are  fairly  severe  (4.2  ft./sec.  RMS  at  400  kts..  6.3  ft./sec.  at  840  kts..  and  8.4  ft./sec.  at  1000  kts.)  and  would  be 
encountered  -nly  occasionally  during  flutter  testing.  In  order  to  give  some  insight  into  the  range  of  answera  which  might  be  obtained 
with  model  matching,  results  are  presented  for  each  of  4  record  segments  (windows)  corresponding  to  each  peak  in  the  response  tracts. 
For  the  two  lower  speeds,  four  windows  were  also  chosen  to  include  the  two  response  peaks  which  occur  in  each  sweep.  Some  of  these 
segments  were  intentionally  chosen  to  be  bad  selections.  For  each  of  the  segments,  frequency  and  damping  values  are  given  for  2nd. 
3rd.  4th,  5th  and  6th  order  models.  The  system  used  to  generate  the  data  is  actually  5th  order  but.  depending  on  the  particular  record 
segment  in  question,  an  adequate  match  may  be  accomplished  with  a  lower  order. 

The  results  given  in  the  table  accompanying  Figure  4  demonstrate  that,  for  a  straightforward  case  such  a»  this,  model  matching 
readily  calculates  an  excellent  damping  and  frequency  value  for  each  of  the  separated  peaks;  slightly  better  answers  arc  obtained  for  3rd 
order  analysis  than  for  2nd  order  Higher  order  models  produce  the  correct  results  for  record  segments  containing  one  peak  but  the 
additional  calculated  roots  arc  difficult  to  interpret.  The  results  for  record  segments  containing  both  peaks  are  naturally  poor  for  2nd 
and  3rd  order  models  but  give  good  results  for  higher  orders. 

The  effect  of  noise  can  be  appraised  by  comparing  the  table  of  Figure  5  with  that  of  Figure  4.  A  slight  degradation  in  accuracy  has 
occurred  and  only  in  the  higher-order  models.  This  noise  case,  an""  the  two  subsequent  ones,  were  analyzed  after  the  data  was  passed 
through  digital  bandpass  filters  with  a  12/db  octave  roll-off  characteristic.  The  corner  frequencies  were  set  at  values  corresponding  to 
the  shaker  instantaneous  frequency  at  the  beginning  and  end  of  each  record  segment. 


• .   ■ 
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Figures  6  anj  7  present  the  results  at  840  kts.  As  above,  3rd  order  results  are  generally  better  than  2nd  order  when  a  single  peak  is 
chosen  for  study  and  the  4th,  Sth  and  6th  order  results  are  generally  better  than  at  400  kts.  Noise  has  degraded  the  data  more  than  at 
400  kts.  It  is  noted  that  the  real  root  listed  in  the  tables  is  seldom  accurate;  an  explanation  for  this  is  that  it  is  necessary  to  choose  a 
window  at  the  beginning  of  the  shaker  sweep  to  obtain  this  root.  Although  no  data  is  presented,  good  results  have  been  obtained  when 
this  is  done. 

The  tables  with  Figures  8  and  9  give  results  at  1000  kts.  for  two  types  of  shaker  forces;  one  which  has  a  slight  variation  in  vector 
amplitude  with  frequency;  and  the  other  with  a  crossover  frequency  essentially  the  same  as  the  coalescence  frequency.  Analysis  of  the 
system  at  this  speed  would  be  a  very  severe  test  for  any  flight  flutter  testing  method,  particularly  with  the  second  force  input,  which 
makes  the  single  peak  appear  to  be  two  separate  peaks. 

Only  4  record  segments  were  chosen  at  this  speed  because  one  true  peak  appears  instead  of  two.  The  tables  show  that  more  scatter 
occurs  in  these  answers  than  at  the  lower  speeds;  the  scatter  is  mainly  in  the  heavily  damped  mode;  the  damping  of  the  critical  mode  is 
fairly  well  predicted.  Other  analyses,  which  are  not  reported  here,  have  also  indicated  that  whenever  two  modes  couple  and  one  is 
lightly  damped  and  the  other  is  heavily  damped,  the  model  matching  analysis  gives  acceptably  good  results  for  the  lightly  damped 
mode. 

Figure  10  gives  results  at  1000  kts  in  the  presence  of  rather  severe  atmospheric  tuibulence.  inspection  of  the  response  trace  shows  a 
considerable  distortion  of  the  signal.  As  expected,  the  results  are  somewhat  inconsistent. 

Some  of  the  record  segments  used  in  the  above  examples  are  not  repre'  ••  i>-  m  of  what  would  be  used  in  actual  practice.  The 
preferred  ones  are  double  starred  at  the  left  of  each  table.  The  answe's  <•-  il  rse  segments  have  been  plotted  in  Figure  3.  The 
resulting  bands  represent  the  degree  of  scatter  which  might  be  expected  using  model  matching.  It  ca  t  be  seen  that  as  the  modes  become 
more  coupled  and  one  gets  closer  to  the  coalescence  speed,  the  scatter  band  becomes  wider.  It  is  the  authors’  view,  however,  that  these 
results  are  at  least  as  accurate  as  those  which  could  be  obtained  by  other  methods  now  in  use. 

3.4  Typical  Flight  Test  Results 

After  gaining  confidence  that  the  model  matching  program  would  predict  accurate  damping  and  frequency  values,  the  program  was 
used  to  do  actual  flight  flutter  testing.  The;  final  evaluation  war  accomplished  by  checking  the  results  with  those  obtained  by  the  old 
shake-and-stop  method.  Satisfactory  agreement  was  obtained. 

Figure  1 1  presents  the  model  matching  program  results  for  a  vertical  tail.  These  results,  together  with  comparable  data  for  the  wing 
and  horizontc  tail,  were  obtained  in  only  six  flights;  this  is  nine  less  than  would  have  been  required  to  obtain  this  quantity  of  data  in 
the  past. 

4.  CONCLUSION 

Although  flight  flutter  testing  methods  in  current  use  in  the  United  States  are  by  no  means  standardized,  they  represent  significant 
advances  over  those  employed  in  the  past,  it  is  rather  encouraging  to  note  that  more  sophisticated  analysis  techniques  producing 
essentially  real  time  results  arc  now  being  widely  used.  The  advances  result  from  the  availability  of  large  high  speed  computers  which  are 
becoming  an  integral  part  of  flight  test  data  management.  Much  more  remains  to  be  done  to  further  improve  present  methods  and  to 
develop  new  ones.  With  a  more  universal. use  of  dedicated  computers,  we  should  expect  a  continuing  proliferation  of  rapid  and  safe 
testing  techniques. 
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Coalescence  Frequency;  No  Gust  Force. 
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Fig.  10  Total  Force  and  Velocity  Response  of  a  Two-Mode  Aeroelastic  Model  at  the  Coalescence  Speed; 
Total  Force  Is  Comprised  of  a  Vane  Force  and  a  Gust  Force  with  an  RMS  Intensity  of  8.4  ft/sec; 
1000  Kts  (94%  Vp)  at  9000  ft;  Vane  Unbalance  Chosen  to  Produce  a  Crossover  Frequency  Near 
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